Intermittency via Self-Similarity — An Analytic Example 
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Turbulence is known to show intermittency. That is, statistical properties vary with the length 
scale in a way not accounted for by statistical similarity where dimensionless ratios of moments are 
constant. Intermittency occurs even in the inertial range of isotropic turbulence, where physical 
intuition calls for a self-similar scale dependence. Perceived as a lack of overall scaling invariance, 
inertial range intermittency has become known as anomalous scaling. We present an analytic exam- 
ple demonstrating how anomalous scaling and self-similarity in the form of global scaling invariance 
can coexist within the same statistics. Whether we observe anomalous scaling or self-similarity 
depends on which variables we consider. Our example illustrates consequences of a symmetry, but 
is not meant as an intermittency model. 

PACS numbers: 47.27.Eq, 47.52. +j, ll.80.Cr, 47.27.Gs, 05.45.Jn, 47.27.Jv, 47.27.eb, 47.27.ed 



Intuitively, fundamental insight is lacking if 'anoma- 
lous' best describes turbulence in its most ideal form. 
Theoretically, turbulence is ideal when all simplifying cir- 
cumstances are present. First, the fluid is incompressible 
and has constant properties-so the motion obeys the in- 
compressible Navier-Stokes equations. Second, there are 
no boundaries and the turbulence is homogeneous and 
isotropic-so there is statistical independence of location 
and orientation. Third, steady large scale forcing puts 
the turbulence in equilibrium and makes its properties 
independent of time. Finally, the Reynolds number is 
high so that a wide range of scales separates forcing at 
the long lengths from dissipation at the short lengths. 
Inertial forces alone govern the physics in this "inertial 
range." Physical intuition dictates that the scale depen- 
dence in the inertial range must be self-similar-in some 
global way. That is, full knowledge of the statistics at 
one scale should suffice to generate the statistics at all 
other scales. 

Kolmogorov 1 (K41) proposed a similarity like that in a 
Brownian motion where dimensionless ratios of moments 
are constant. Applied to the difference 5v(£) in a veloc- 
ity component over a distance £, this so-called statistical 

similarity requires (Sv p ) / (8v 2 ^ p ^ 2 to be independent of 
I as suggested by dimensional analysis. Statistical simi- 
larity, however, does not agree with the evidenced The 
moment ratios are not constant. Equivalently, the prob- 
ability density functions (pdf's) for Sv(£) at various I do 
not collapse to one when plotted in units of standard 
deviations. Instead, one observes intermittency: fluctu- 
ations of many standard deviations become increasingly 
more likely at smaller I. Correspondingly, the tails of 
the pdfs flare out and become thicker. Although the mo- 
ments are power laws, e.g., (Sv p ) oc !!•* , the exponent 
£ p varies non-linearly with the order p. In contrast, sta- 
tistical similarity requires that £ p be linear. The name 
"anomalous scaling" alludes to historical frustrations in 
understanding why statistical similarity fails. Frisch 2 
presents an account of the history and theories for in- 
termittency. Kolmogorov 3 addressed the issue with his 
log- normal model (K62). Although K62 has fallen out 



favor, it remains among the best-known intermittency 
models. To a large extent, it accounts for the discrepan- 
cies between the K41 predictions and experimental ob- 
servations; see Ref.0, p. 258. However, new models, e.g. 
Refs.|5j|a, describe observations better than K62. 

Here, we show how self-similarity and anomalous scal- 
ing can coexist within the same statistics. The variables 
and functions we choose determine whether we observe 
self-similarity or only anomalous scaling. In particu- 
lar, self-similarity emerges from considering the energy 
at scale I. To avoid lengthy mathematical derivations, 
we illustrate these ideas by an example. The example 
comes from a similarity theory we have developed for 
the inertial ranged This theory puts the inertial range 
pdf for Navier-Stokes turbulence within the small class 
of functions generated by the inverse Mellin transform 
of z^ 1 exp (sign (/? — 1) z 13 ). One of these functions has 
K62 scaling exponents, namely (3 — 2. We choose this 
function in our example for three reasons. First, the cal- 
culations work analytically and in a closed form indepen- 
dent of the theory J Second, K62 is widely known and is 
in reasonable agreement with observations. Third, the 
concepts we emphasize are shared by all functions in the 
class: global scaling invariance (self-similarity) emerges 
through the energy at scale £ and coexists with intermit- 
tency within the same statistics. 

Our example uses Mellin transforms extensively, so we 
review the needed properties. The Mellin transform 
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is useful for dealing with moments on the positive real 
axis. If known, $(z) provides all moments including frac- 
tional orders. On the full axis, we have 



x p <p(x)dx = M [<j>(x);p + l] + (-l) p 371 [<f>(-x);p + 1] . 

. (2) 

However, full range moments are defined only for inte- 
ger orders. Consequently, we face "Hausdorff's moment 
problem," where the integer moments do not uniquely 
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identify a function; see Ref.|8|, p. 21. For example, many 
functions have the same integer moments as the log- 
normal. In contrast, the inverse Mellin transform is 
unique. The Mellin transform is closely related to the 
Fourier transform; see Ref. 9, p. 41. Each operational rule 
for the Fourier transform has a counterpart in the Mellin 
transform. We need the ruleitt 
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(j){{x/a) 1/q ^ ;z = a z q${qz), a,q>0. (3) 



From tables of transforms^ we need the entry 
9Jt [| eric In r) ; z] = z^e*' , 



(4) 



where erfc(t) = 2tt~ 1 / 2 J" t °° e~ T dr is the complementary 
error function. 

For our example, we design a function F(r,£) with 
global scaling invariance in p = lnr, but with anoma- 
lous scaling in r. The self-similarity F(r,£) = f({p — 
p(£))/a(£)) expresses our scaling invariance as an afhnc 
transformation on the p— axis. That is, p and a reset 
the origin and the unit on the axis in the same way as 
we traditionally use mean and standard deviation. We 
shall call the invariance under this affine transformation 
'normal scaling' to draw the distinction from 'anomalous 
scaling.' Anomalous scaling requires the p'th moment of 
F(r, £) to be a power law C p £^ p where £ p depends non- 
linearly on p. Only specific combinations of f(p), p{£), 
and a(£) permit F(r, £) to have anomalous scaling in r. 
f(p) = i erfc (^p) is one such function as we shall now 
show. By construction, we have 



1/<T 



(5) 



Using J2J we compute the moments: 

/ r p F(r,£)dr = Wl[F(r,£);p + 1] 
Jo 

= ae {p+1)fi M [/(lnr); (p + 1) a] = 



p+1 



(6) 



For £ < 1 (small scales), we choose p(£) — — ahxt and 
a(£) = jyj— ln£ so as to obtain the power laws 



tf-(p+l)a-0+l)V 

r p F(r,£)dr = . (7) 

o P+ 1 



So, the scaling in r is anomalous except when 7 = 0. 

Let us establish a connection with turbulence. Be- 
cause of incompressibility, the velocity is divergence-free 
and consequently has only two independent components. 
Thus, two random variables describe isotropic turbulence 
at each scale. One choice is longitudinal and transverse 
velocity increments, i.e., Svu{£) and Sv±(£). Another 
is left- and right-handed velocity components obtained 




FIG. 1: (Color online) The joint pdf J'(x,y,t) is shown for 
I = 0.01 together with the one dimensional pdf X(x,£) and 
y(y,£). The radial profile P(r,l) is self-similar (Fig. 2), but 
X(x,l) shows anomalous scaling (Fig. 3). 



via the complex helical waves decompositioniii Theoret- 
ically, we want the left-right symmetry of the Navier- 
Stokes equations reflected in the variables. Thus, their 
squares summed with equal weight should yield the en- 
ergy in scale £. That then provides a Cartesian descrip- 
tion. To emphasize this point, we call the two random 
variables X and Y. In this regard, left- and right-handed 
amplitudes are good variables?^ but Sv\\(£) and Sv±(l) 
are not; see Ref.ll2l p. 208. Let X(x,£) be the pdf for 
X, i.e., X{x)dx = Pr{x < X < x + dx}. Similarly, let 
y(y, £) be the pdf for Y, and J{x, y, £) the joint pdf for X 
and Y. Because J{x, y, £) — > as x 2 +y 2 — > 00, the graph 
of J has the shape of a mountain centered near (0,0). 
Thus, we use polar coordinates (x,y) — r (cos 0, sin 0) 
and we obtain an azimuthal decomposition 

J(x,y) = P(r)+cos{0-0 1 )P 1 {r)+cos{2(0-0 2 ))P 2 (r) + - ■ ■ 

(8) 

where 61,62, ... are phase-constants. We focus on the ax- 
isymmetric contribution to J . Returning to our example, 
we construct J from F(r, £) : 



J(x,y,£) = P(r,£) = C(£)F(r,£), 



(9) 



where the factor C{£) = 1/ (2tt971 [F(r, £); 2]) = 
7r~ 1 e~ 2M-4<T gives unit volume under the graph and en- 
sures J is a pdf. 

With normal scaling built into J , consider the situa- 
tion graphically. Figure 1 shows J{x,y,l) for £ = 0.01. 
As we vary £ the graph changes. p(£) controls the half- 
width-radius of the "pie," while a(£) sets the slope there. 
p(£) and a(£) do not change proportionally, so statistical 
similarity is impossible. In the limit £ — > 1 — , the graph 
looks like a disk, but as £ — > + a sharp peak develops at 
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FIG. 2: (Color online) (a) The normalized radial profile 
P(r,£)/C(£) = F{r,£) for i = 0.99, 0.5, 0.1, 0.01, 10 -20 ; 
(b) Same graphs on a logarithmic abscissa p — lnr; (c) Col- 



lapse by normal scaling p - 
heavy black graph). 



(p-ix(i))/a(£) to ierfc(ip) (the 



(0,0). Figure 2a shows the normalized radial profile for 
various I. For small £, the curves cluster near the verti- 
cal axis and are impossible to distinguish. The remedy 
is a logarithmic axis (Fig. 2b). On the logarithmic axis, 
the curves look alike so we align the midpoints (Fig. 2c). 
Then horizontal scaling (i.e. multiplying the abscissa by 
1/ct) collapses all curves onto f(p) = \ erfc i\p) ■ 

In contrast, the traditional analysis of the same statis- 
tics shows anomalous scaling. Two steps are involved. 
First, we find scaling exponents £ p for the structure func- 
tions (i.e., moments) of a single random variable such as 
X. Second, we plot the corresponding pdf, X(x,£), for 

1/2 

various £ with the abscissa in units of {X 2 ) and the 
ordinate scaled to give unit area. 

To obtain £ p , we calculate the moments of X(x,£) by 
integrating J(x,y,£) over y. Using (JJJ and © we have 



/>oo />oo 

x p X(x,£)dx = / x p / P(r,£)dydx 

II JO J-QO 

7r/2 />oo 

/ (r cos d) p P{r,£)rdrd9 

-7T/2 J-OO 

/>7r/2 />oo 

= 2/ cos p ddd / r p+1 P(r,£)dr 
Jo Jo 

K p Tt[P(r,£);p + 2} = K p r ap - ( - p2+4p ^ 2 (10) 



where p > -1 and K p = (§ + \) /V (f + l) with 
T(z) being the usual gamma function. (Full range mo- 
ments follow from Note the 9— integral in i|10|) 
converges only when p > — 1, whereas the r integral 
converges for p > —2. Thus, by considering X rather 
than \/X 2 + Y 2 we loose all exponents in the range 
-2 < p < -1. Clearly, C P = -ap - {p 2 + ^ph 2 - So 
Co = and C3 = -3a - 2I7 2 . Choosing a = -§ - 77 2 
and 7 = 1/V90 yields ( 3 = 1 and Ce = 1.8. By de- 
sign, our £p, shown in the inset in Fig. 3, have classical 
K62-valuesi£ 

Figure 3 shows the graphs of X(x, £) for the same val- 
ues of £ as in Fig. 2. This time, the curves do not collapse 
to one, but exhibit the usual flaring tails associated with 
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FIG. 3: (Color online) (Inset) K62 Scaling exponents; (Main 
figure) X(x,£) obtained numerically from 112H for the same 
values of I as in Fig. 2. 



anomalous scaling. To reveal asymptotic properties, we 
rewrite X(x,£) analytically. For x > 0, 

/oo poo 
J(x,y,£)dy= P(r,£)dy 
-OO J — OO 

fir II 

= / P (x sec 9, 1) x sec 2 9 d9 

J -it/2 

2x c*/ 2 r°° 2 
= -=C(£) I I e _< sec 2 9 dtd9 



xC 



ir/2 



exp 



^ ln(e — sec 9) 

In 2 (e _Al xsec( 
i^ 2 



tan 2 9d9, (11) 



where the last step requires integration by parts. Upon 
the substitution rj = Insect, the integral becomes nu- 
merically friendly with nice asymptotic properties: 



X(x,£) 



c r V 7r Jo 



v/l-e-^exp^-^^j ) dr > 

where r\ c = 2a 2 —Inx + fx. For fixed £, we have rj c — > oo 
as x — > + and the integrand reduces to a Gaussian so 
that 

#(G, £) = ^ V4tt(t 2 = 27r- 1 e- 3CT2 ^. (13) 

With (X 2 ) 1 ' 2 = ±e 6CT2 +^, we have (X 2 f 2 X(Q, £) = 
^-l^ 2 = ^-1^-1/30 _> oo as £ ^ 0+. Thus, the peak 
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In (K(r,^)/r) 




FIG. 4: (Color online) (a) The pdf for the amplitude 
V 'X 2 + Y 2 at scale £ with the same £ values as in Figs. 2 and 
3; (b) lZ(r,l)/r — 2nP(r,t) is self-similar and collapses by 
normal scaling; (c) Structure functions S p (£) = 2£^ p /(p + 2) 
plotted to reveal the virtual origin for the inertial range. 



in Fig. 3 rises without bound as £ — > + . Again for fixed 
£, we have r\ c — ► — oo as x — > oo so that the essential con- 
tribution to the integral (|12[) comes from a small interval 
to the right of x = 0. Asymptotic analysis yields: 



*~^ e - 4 ^exp 



ln 2 (xe-^) 



7T " \ 4CT 2 



1 -3/2/ 

In ' (a 



2 7 2 |ln 



as x — > oo. The decay of (cc, ^) is essentially log-normal. 

We emphasize that our example is not intended as 
a model of intcrmittency, but rather as an illustra- 
tion of the symmetry expressed by the global invariancc 
of F(r,£) = P(r,£)/C(£) under the affine transforma- 
tion p(r) = lnr — > (p — p(£)) / a{£). Our example is 
not a model because we make no modeling approxima- 
tions. Rather, we suggest that equilibrium turbulence 
has precisely this symmetry. Considering Navier-Stokes 
equations with steady forcing as a dynamical system, 
equilibrium turbulence represents an attractor in phase 
space. Our suggestion is that Navier-Stokes equations 
restricted to the living space of the attractor obey our 
symmetry. Our suggestion is supported by shell model 
calculations^ and recently by low resolution DNS in a 



periodic box. Theoretically, we can show 7 that the sym- 
metry exists together with anomalous scaling in r only 
when 37t[/(lnr); z] — z" 1 exp(sign(/3 — l)z@). Our pre- 
liminary DNS show « 1.4, significantly less that two. 
Thus, there is no point in attempting to fix and resurrect 
K62. 

Using our example, we illustrate how our global scal- 
ing invariance appears in analysis of data with J(x, y, £) 
statistics. With X 2 + Y 2 being a random variable for the 
energy at scale I, we construct the pdf for the velocity 
magnitude as 7l(r, £) dr = Pr{r < V X 2 + Y 2 < r + dr}. 
A simple and robust technique for that purpose is kernel 
density estimation^ Plotting the pdf on log-log scales 
for various £, we obtain Fig. 4(a). All curves have an 
asymptote with slope one. This is the signature of a 
monopole in polar coordinates, i.e. < P(0, £) < oo 
implies TZ(r,£)dr w 27rP(0, £) r dr as r — > + . We make 
the asymptote horizontal and thus immune to abscissa 
scaling by considering the function 1Z(r,£)/r, which is 
2nP(r,£) so that the curves in Fig. 4b collapse to one 
by normal scaling. In addition, our global scaling invari- 
ance also shows up in S p {£) = ((V X 2 + Y 2 ) p ). Using (JSJ) 
and ©, we have S p {£) = 2f"/\p + 2). When plotted as 
in Fig. 4c, all lines converge on a focal point. It is at 
(0, 0) because of the way we have normalized our exam- 
ple. The focus marks an infrared virtual origin for the 
inertial range. The focal abscissa is the intrinsic length 
scale for inertial range. The ordinate is proportional to 
the dissipation as in Kolmogorov's four-fifths law, but we 
do not yet know the proportionality constant. 

Consider the juxtaposition of Figs. 2 and 3, which 
represent the same statistics, J{x,y,£). In Fig. 2 all 
curves collapse to one by normal scaling, while in Fig. 
3 anomalous scaling prevents them from doing so. Our 
example is analytic and exact-no modeling or numerical 
issues blur the juxtaposition. Our example shows directly 
that intermittency is not at odds with global scaling in- 
variance. Normal scaling on the logarithmic axis accom- 
modates two independent parameters p (£) and a {£) ; see 
Fig. 2(b)-(c). In contrast, statistical similarity allows 
only one and fails for that reason. 
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